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Quadrupolar Frustration in S hast ry- Sutherland Lattice of DyB4 Studied by Resonant 
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We have observed geometrical frustration of quadrupolar and magnetic moments in dyspro- 
sium tetraboride, DyB4, where the rare-earth sites form a Shastry- Sutherland lattice. Resonant 
X-ray scattering at the L\\\ absorption edge of Dy was utilized. Analysis of the energy, polariza- 
tion, temperature, and azimuthal-angle dependences of the El resonance of the (10 0) forbidden 
reflection show that the magnetic and quadrupolar components within the frustrated c plane 
have a short-range correlation, suggesting that the moments are fluctuating. In contrast, the 
basic antiferromagnetic component along the c-axis has a long-range order. 

KEYWORDS: geometrical frustration, DyB4, resonant X-ray scattering, quadrupole order, fluctuation 



Magnetic materials with geometrically frustrated mag- 
netic interactions exhibit a wide variety of intriguing phe- 
nomena caused by fluctuating spins that cannot order 
down to very low temperatures. For example, hexagonal 
ABX3 (A=alkali metal, B=transition metal, X=halogen) 
compounds such as CsCoCla with a triangular lattice of 
Ising spins exhibit an unusual antiferromagnetic state in 
which one of its three sublattices is disordered.^ Frustra- 
tions in kagome and pyrochlore lattices offer a platform 
where a macroscopic number of spin configurations are 
degenerate, giving rise to short-range correlated fluctuat- 
ing spins called spin liquid states.^ These exotic phenom- 
ena are caused, of course, by the spin degree of freedom 
of the magnetic ions. On the other hand, there are only 
a few examples that suggest frustration of the orbital 
degree of freedom. In LiNi02 and FeSc2S4, for example, 
their anomalous magnetic properties are interpreted in 
association with the orbital frustration.^'^ However, to 
the best of our knowledge, there has been no experiment 
in which the frustrated orbitals were observed directly 
by an orbital-sensitive method. In this letter, we present 
direct evidence for fluctuating quadrupolar and magnetic 
moments in a rare-earth compound, DyB4, which has a 
Shastry- Sutherland- type geometrical frustration.^ 

Rare-earth tetraborides, RB4, have been attracting 
growing interest as a system in which the rare-earth net- 
work involves a geometrical frustration. The R ions are 
at the sites of the tetragonal space group of PA/mhm: 
(^^^^1 0), R2 (|-x,x,0), |-x,0), 
and i^4 = (I + X, -x,0), with x = 0.3175 for DyB4.^ 
This is equivalent to the Shastry- Sutherland lattice as 
illustrated in the inset of Fig. l(b).^ It is noted that 
the nearest-neighbor and next nearest-neighbor distances 
differ only by 0.33% and that the Dy^+ ion has a 
huge magnetic moment of lO/iB- Therefore, the Shastry- 
Sutherland lattice of Dy can be regarded as a combi- 
nation of triangles and squares of classical moments. 
In addition, an orbital degree of freedom is active in 
DyB4 as described next. Consequently, the way of remov- 
ing frustration in DyB4 is expected to be different from 
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the orthogonal dimer system of quantum spins such as 
SrCu2 (303)2, where the frustration can be removed by 
making spin-singlet dimers.^'^ 

The physical properties of DyB4 have been studied in- 
tensively by Watanuki et al.^'^ Two phase transitions 
take place, at Tni = 20.3 K and at Tn2 = 12.7 K. Below 
Tni, it is established by neutron powder diffraction that 
an antiferromagnetic order takes place with the magnetic 
moments aligned along the c axis; those on Dy(l) at 
Ri and Dy(2) at R2 are directed upwards and those on 
Dy(3) at R^ and Dy(4) at R/^ downwards. As a result, 
forbidden (odd 0) reflections arise, but unit cells do 
not change. What is intriguing is that the C44 mode of 
the elastic constant keeps softening even below Tni- The 
softening, which follows the normal Curie law in phase 
I (T > Tni) that indicates a quadrupolar degeneracy, is 
even more enhanced in phase II (Tn2 ^ T < Tni). More- 
over, in phase II, strong ultrasonic attenuation occurs, 
suggesting fluctuation of the quadrupolar moments. Spe- 
ciflc heat measurement shows that the entropy of In 2 
and i?ln4 is released at Tn2 and at Tni, respectively, 
with increasing temperature. Therefore, the ground state 
is a pseudo-quartet with quadrupolar degeneracy, which 
is somehow not lifted even in phase II. In phase III 
(T < Tn2), the ultrasonic attenuation and softening stop; 
neutron powder diffraction suggests that a c-plane mag- 
netic component arises in phase III.^^ These anomalous 
behaviors strongly suggest that phase II is not a normal 
antiferromagnetic state and that there remains fluctuat- 
ing quadrupolar and magnetic moments, as flrst pointed 
out by Watanuki et al7 

In order to investigate phases II and III of DyB4, we 
have utilized resonant X-ray scattering (RXS), which 
has been developed as a powerful tool to probe both 
quadrupolar and magnetic moments. Since RXS mea- 
sures the average structure of the relevant moments in a 
very short time scale determined by the lifetime of the 
intermediate state (~ 10~^^ sec), it would be possible 
to observe a snapshot of the correlated moments if the 
fluctuation rate is slower than the time scale of the ob- 
servation. In the present experiment, we have detected a 
broadened forbidden reflection, indicating a short-ranged 
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Fig. 1. Incident energy dependences of the peak-top intensity of 
the (1 0) reflection for ct-tt' at (a) ^ = 0° and (b) ^ = 90°. 
The peak around 7.824 keV for ^ = 0° is due to the multiple 
scattering. The top figure shows the fiuorescence spectrum. The 
inset shows the Dy lattice; the nearest neighbor bonds are shown 
by the thick lines. 



quadrupolar and magnetic order. 

A single crystal of DyB4 was grown by the floating 
zone method using a high-frequency furnace. The sample 
quality was checked in terms of magnetic susceptibility 
and electrical resistivity, which showed good agreements 
with the results in ref. 6,7. Synchrotron X-ray scattering 
experiments were performed using a four-circle diffrac- 
tometer at beamline BL-16A2 of the Photon Factory in 
KEK. The incident X-ray was tuned near the Lni ab- 
sorption edge of Dy (~ 7.79 keV), where the resonance 
of 2p ^ 5d dipole transition occurs. We investigated the 
{h 0) reflections using a sample with the (100) surface. 
Azimuthal-angle scans and polarization analysis using a 
PG analyzer-crystal were also carried out; the azimuthal 
angle ^ is deflned to be zero when the c-axis lies in 
the scattering plane. The mosaic width of the sample 
was about 0.03° full width at half-maximum (FWHM) 
in phase I, indicating the high quality of the crystal. 

Figure 1 shows the incident energy dependences of the 
(10 0) forbidden reflection for the ct-tt' channel at 6 K 
(phase III), 15 K (phase II), and 25 K (phase I) at two 
azimuthal angles of ^ = 0° and 90°.^^ The resonant 
enhancement is clearly observed around the absorption 
edge. We notice three resonant features at ^=0°, which 
are indicated by an arrow at 7.79 keV and two bars at 
7.795 keV and 7.801 keV. The resonance at 7.79 keV, 
the inflection point of the fluorescence spectrum, disap- 
pears above Tni; this can be ascribed to the resonant 
exchange magnetic scattering originating from the 4/ 
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Fig. 2. Azimuthal-angle dependence of the integrated intensity 
of the (1 0) resonant reflection at E = 7.79 keV for a-ir' . The 
lines are fits with cos-^ ^ and sin^ ^. 



magnetic moment. The other two resonances at 7.795 
keV and 7.801 keV, which exist even above Tni, can be 
ascribed to the anisotropic tensor susceptibility (ATS) 
scattering that reflects the local crystal-fleld anisotropy 
of the surrounding boron atoms. We call this boron- ATS 
scattering hereafter. The double peak energy-spectrum 
and the cos^ ^ dependence in intensity are the same as 
those in GdB^.^^' At ^=90°, the boron- ATS scattering 
vanishes and only the resonance at 7.79 keV is observed 
as indicated by the arrow. We also observed nonresonant 
reflection below the edge, which we ascribe to magnetic 
scattering; the intensity is proportional to the resonance 
at 7.79 keV. 

We note that the most mysterious point is the ex- 
istence of the intensity at ^=90° in phase II. Accord- 
ing to the simple magnetic structure deduced from neu- 
tron diffraction, where the antiferromagnetic moments 
align completely along the c axis,^ the structure factor of 
the resonant magnetic scattering becomes proportional 
to cos^. However, this is inconsistent with the data of 
Fig. 1(b) for 15 K. The azimuthal-angle dependence of 
the resonance intensity shown in Fig. 2 also demonstrates 
that the intensity at ^=90° does not vanish in phase 11. 
In phase III, a complex magnetic structure is proposed 
by neutron scattering, where the c-plane component of 
the magnetic moments appears.^ It is therefore expected 
that the c-plane component gives the intesity at ^=90°. 

Figure 3 shows the temperature dependences of the 
integrated intensity and half width at half maximum 
(HWHM) of the (1 0) resonant reflection at 7.79 keV. 
There is an obvious difference between the behavior for 
ilr = 0° and 90°. At ^ = 0° the intensity rises up steeply 
below Tni and saturates below about 17 K, whereas at 
^ = 90° it increases gradually below Tni and exhibits 
a sudden increase below Tn2- Moreover, the HWHM at 
^ = 90° is obviously broad and exhibits a sudden drop 
below Tn2, whereas at ^ = 0° it is almost resolution 
limited and exhibits a slight broadening in phase 11. It 
is noted that the rocking (u) scans in Fig. 3 correspond 
to the /-scan at ^ = 0° and /c-scan at ^ = 90° in the 
hkl reciprocal space. On the other hand, the HWHM for 
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Fig. 3. Temperature dependences of (a) the integrated intensity 
and (b) the HWHM of the rocking scan (uu) of the (1 0) res- 
onant reflection at ^ = 0° and 90° {E = 7.79 keV, a-n'). The 
solid lines are guides for the eye. The dotted line indicates the 
resolution limit estimated from the fundamental reflections. The 
inset shows the peak proflle at 7 K and 15 K for ^ = 90°. 



the /-scan at ^ = 90° (-0.01 A"^ at 16 K and -0.007 
at 6 K) is also obviously broader than the resolution 
(-0.005 A~^), whereas that for the A:-scan at ^ = 0° is 
resolution limited, although these resolutions are much 
worse because the scans correspond to the x-scan. These 
results suggest that a different order parameter exists at 
^=90°, which has a short correlation length in phase IT 

We briefly refer to the structural phase transition in 
phase III that has been studied using an X-ray from a 
Mo target (17.48 keV). The single peak of the (0 6) re- 
flection at 16 K splits into four peaks along the h and k 
directions at 9 K; the (10 0) peak becomes broad along 
the / direction. These results, together with the survey of 
some other Bragg peaks, show that the structure changes 
from tetragonal to monoclinic. The angle Zac, or Z6c, is 
estimated to be 89.84° at 9 K. It should be remarked 
that the antiferromagnetic domains with the (10 0) and 
(0 1 0) propagation vectors result in the monoclinic dis- 
tortions within the be plane {Zbc ^ 90°) and ac plane 
(Zac 7^ 90°), respectively; this is inferred from the fact 
that the (10 2) resonant magnetic Bragg peak splits only 
along the k direction in phase III. 

Let us discuss the origin of the nonzero intensity of 
the broad (10 0) resonant reflection at ^ = 90° in phase 
II. We analyze the RXS result by taking into account 
the c-plane components of both magnetic (rank 1) and 
quadrupolar (rank 2) moments. For this purpose, we use 
a theory developed by Lovesey et al. which directly con- 
nects the atomic tensors to the scattering amplitude of 

Here, we examine two possible models of the in-plane 
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Fig. 4. Models of possible magnetic and quadrupolar structures. 
The local x- and 2:-axes correspond to the c-plane and c-axis 
magnetic components of each Dy ion, respectively. Ellipses with 
light and dark shades represent {Ozx)-type charge distributions 
that are extended above and below the plane of the paper, re- 
spectively. 



magnetic and quadrupolar structures that are consistent 
with the experimental results. They are illustrated in 
Fig. 4. We take the directions of the c-plane and c-axis 
components as the local x- and z-axes, respectively. The 
atomic tensors expressed in the local xyz coordinates are 
assumed to be the same for all the four Dy ions; they can 
be transformed to each other by appropriate Euler rota- 
tions. Quadrupolar moments represented by the ellipses 
can coexist because of the strong spin-orbit interaction. 

The unit cell structure factor of the {h=odd 0) re- 
flection for the two models in Fig. 4 is written as 

F^.^V = 2A:i6cos6>(V2(J^)cos^-(J^)sin^) 

T2k2a cos i9( ^2(022) cos ^ + {O^x) sin ^) (1) 

for El resonance and the ct-tt' channel, where the minus 
sign in the second line corresponds to model I and the 
plus sign to model II; a = cos(27r/ix) and b = sm{27rhx); 
is the Bragg angle; ki and k2 are constant real num- 
bers. Note that the tensors in eq. (1) are of the 5d state; 
(Ja) can be induced through magnetic exchange interac- 
tion with {Ja)Af] through d-f Coulomb interaction 
with {0^)4f or through mixing with the 2p orbitals of 
boron atoms (boron- ATS mechanism). The boron- ATS 
mechanism obviously gives rise to (O22), and this term 
explains the azimuthal- angle dependences of the reso- 
nances at 7.795 keV and 7.801 keV in phase I. In phase 
II, the c-axis component of the 4/ magnetic moment in- 
duces (Jz). The {Jz) term explains the cos^ ^ component 
in the azimuthal-angle dependence of the resonance at 
7.79 keV as shown in Fig. 2. 

The {Jx) and (Ozx) terms, with the sin^ dependence, 
can explain the nonzero intensity at ^=90°. However, 
because the cos ^ and sin ^ terms interfere coherently in 
eq. (1), the structure factor can be transformed into the 
form Acos(^ + a), where A and a are constants. This 
means that the intensity must vanish somewhere in the 
azimuthal scan, which is the case in GdB4.^^ However, 
the experimental data in Fig. 2 do not support this sim- 
ple scenario. The data suggest that the cos^ and sin^ 
terms contribute independently rather than interfere co- 
herently. This is demonstrated by the lines in Fig. 2. It 
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seems that the resonance at ^ = 90° represents the in- 
plane components {Jx) and (Ozx)^ and the resonance at 
^ = 0° represents the (J^) and (O22) moments, respec- 
tively. 

A possible reason for this behavior is that the cor- 
relation length the of (Jx) and (Ozx) components is 
much shorter than that of (J^) and {O22)' By fitting 
the /c-scan profile at ^ = 90° with a Lorentzian that 
is convoluted with the resolution function, the correla- 
tion length of {Jx) and (Ozx) along the 6-axis is esti- 
mated to be (1.4 zb 0.1) X 10^ A at 15 K in phase II and 
(4.2 zb 0.2) X 10^ A at 7 K in phase III. The correlation 
length along the c-axis is also short, less than the res- 
olution of 1.5 X 10^ A although precise estimation was 
difficult. On the other hand, the correlation of (Jz) and 
(O22) along the c-axis estimated from the /-scan profile 
at ^ = 0° is (9.2 zb 0.5) x 10^ A at 15 K and over 2 x 10"^ 
A, the resolution limit, at 7 K. The correlation along the 
6-axis is also expected to be longer than the resolution 
limit of 1.5 X 10^ A. Then, the volume fraction from which 
the scattering occurs coherently is much smaller for (Jx) 
and (Ozx) than for (Jz) and {O22)' Therefore, the cos^ 
and sin^ terms in eq. (1) almost do not interfere and 
they give reflections independently. 

Thus, we consider that the resonance at ^ = 90° re- 
flects {Jx) or (Ozx)^ which is induced by the magnetic or 
quadrupolar moments of the 4/ electrons, respectively. 
The peak profiles are clearly broadened in phase II, indi- 
cating that the correlation length is short. Furthermore, 
these in-plane components are considered to be fluctuat- 
ing. This is supported by neutron powder diffraction in 
phase II, which shows that {gJz)4f is only 6.9/iB while 
the full moment is 10/iB-^ The strong ultrasonic atten- 
uation in the C44 mode also supports the ffuctuation of 
(Ozx)-^ In phase III, neutron diffraction shows that an 
in-plane component of 5.3/xb appears, resulting in a total 
moment of 8.8/iB; the ultrasonic attenuation and soften- 
ing also stop. RXS at ^=90° shows that the correlation 
length becomes longer, although it is not as long as the 
resolution limit. These results show that the fluctuation 
stops in phase III and a static short-range order of {Jx) 
and {Ozx) remains. 

Which of {Jx) and {Ozx) is dominant in this reflec- 
tion? Unfortunately, it is not distinguishable from the 
present experimental results and from eq. (1) because 
the two factors are completely in phase. In the present 
case of DyB4, because {Jz) is flnite, {Ozx) will be induced 
if {Jx) arises because of the strong spin-orbit coupling, 
and {Jx) will be induced if {Ozx) arises as well. Ultra- 
sonic attenuation and softening of the C44 mode seem to 
support the dominance of {Ozx)- It would be possible to 
distinguish if we could investigate by RXS the relative 
change in the signal for h =1, 3, and 5 of the {h 0) 
reflection through changes in the factors a and which 
have different h dependences. 

With respect to the origin of the short-range correla- 
tion and fluctuation of the in-plane moment, one plausi- 
ble scenario is that the {Ozx) moments in the triangular 
connection of the Dy lattice experience frustration. If we 
assume the charge distribution of the 4/ electrons prefer 
to avoid each other, the arrangement of {Ozx) in models 



I and II experiences frustration, which could lead to an 
orbital liquid state. 

We consider that the static order of the in-plane mo- 
ments takes place when the frustration is removed by 
the lattice distortion. The arrangement of the {Ozx)-^yp^ 
quadrupolar moment in models I and II will favor a uni- 
form lattice distortion through a cooperative Jahn- Teller 
effect into a monoclinic structure where the be plane is 
distorted {Zbc 7^ 90°). This is consistent with the ex- 
perimental result in phase III. A possible space group of 
the monoclinic phase allowed by Landau theory is P2i/c, 
where all the atomic sites are described by the 4e site.-^'' 
It is intriguing that the peak width at ^ = 90° is still 
broader than the resolution limit, which means that the 
correlation length does not diverge. 

In summary, we have investigated the phase transi- 
tions occuring in the Shastry- Sutherland lattice of DyB4 
by resonant X-ray scattering, and detected the short- 
range correlated in-plane moments {Jx) and {Ozx)- They 
are expected to be fluctuating in the high-temperature 
antiferromagnetic phase because of frustration and they 
change into possessing a static order with short-range 
correlation in the low-temperature phase with monoclinic 
distortion. The present results indicate that the doublet 
ground state remains well below T^i despite the mag- 
netic ordering, at which the c-axis component is flxed 
and the pseudo- quartet should be split into four singlets 
by the Zeeman effect. Why does the doublet exist in 
Phase II? What kind of quadrupolar degeneracy is real- 
ized? Such fundamental issues need to be solved in future 
studies. 
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